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ABSTRACT 

We present line profiles and profile parameters for the narrow-line regions (NLRs) of six Seyfert 1 
galaxies with high-ionization lines: MCG 8-11-11, Mrk 79, Mrk 704, Mrk 841, NGC 4151, and NGC 
5548. The sample was chosen primarily with the goal of obtaining high-quality [Fe vn] 26087 and, when 
possible, [Fe x] 26374 profiles to determine if these lines are more likely formed in a physically distinct 
“coronal line region” or are formed throughout the NLR along with lines of lower critical density (n cr ) 
and/or ionization potential (IP). 

We discuss correlations of velocity shift and width with n cr and IP. In some objects, lines of high IP 
and/or n cr are systematically broader than those of low IP/n cr . Of particular interest, however, are 
objects that show no correlations of line width with either IP or n cr . In these objects, lines of high and 
low IP/n cr are remarkably similar, which is difficult to reconcile with the classical picture of the NLR, in 
which lines of high and low IP/n cr are formed in physically distinct regions. We argue for similar spatial 
extents for the flux in lines with similar profiles. 

Here, as well as in a modeling-oriented companion paper (Paper II), we develop further an idea sug- 
gested by Moore & Cohen that objects that do and do not show line width correlations with IP/n cr can 
both be explained in terms of a single NLR model with only a small difference in the cloud column 
density distinguishing the two types of object. Overall, our objects do not show correlations between the 
full width at half-maximum (FWHM) and IP and/or n cr . The width must be defined by a parameter that 
is sensitive to extended profile wings in order for the correlations to result. In Paper II, we present 
models in which FWHM correlations with IP and/or n CT result only after simulating the lower spectral 
resolution used in previous observational studies. The models that simulate the higher spectral resolution 
of our observational study produce line width correlations only if the width is defined by a parameter 
that is more sensitive to extended profile wings than is the FWHM. Our sample of six objects is in effect 
augmented by incorporating the larger sample (16 objects) of Veilleux into some of our discussion. 

This paper focuses on new interpretations of NLR emission-line spectra and line profiles that stem 
directly from the observations. Paper II focuses on modeling and complements this paper by illustrating 
explicitly the effects that spatial variations in electron density, ionization parameter, and column density 
have on model profiles. By comparing model profiles with the observed profiles presented here, as well 
as with those presented by Veilleux, Paper II yields insight into how the electron density, ionization 
parameter, and column density likely vary throughout the NLR. 

Subject headings: galaxies: nuclei — galaxies: Seyfert — line: profiles 


1. INTRODUCTION 

One well-known observational feature of the narrow-line 
region (NLR) in Seyfert galaxies is the preponderance of 
blueward asymmetries in the emission-line profiles (e.g., 
Heckman et al. 1981; Whittle 1985a; Dahari & De Robertis 
1988; De Robertis and Shaw 1990; Veilleux 1991a). To 
account for the asymmetry, it is accepted in general that the 
NLR velocity field is mostly radial and that some source of 
opacity, most likely dust, is present. However, a consensus 
regarding the direction of cloud flow has yet to emerge 
because the blueward asymmetry can be explained qualit- 
atively by either internal obscuration in the case of inflow or 
by external obscuration in the case of outflow (e.g., 
Capriotti, Foltz, & Byard 1981). 

Another well-known observational feature of the NLR is 
that the full width at half-maximum (FWHM) of emission 
lines often correlate with ionization potential (IP) and criti- 
cal density (n cr ). This feature has been commonly used to 
argue for a velocity field increasing radially inward under 

1 Center for Astrophysics and Space Sciences 0111, 9500 Gilman Drive, 
University of California, San Diego, La Jolla, CA 92093-01 1 1. 

2 Department of Astronomy, University of California, Berkeley, Berke- 
ley, CA 94720. 


the assumption that lines of low IP/n cr are emitted less 
efficiently (for fixed emitting mass) close to the central 
source than lines of high IP/n cr (e.g., Osterbrock 1981; De 
Robertis & Osterbrock 1984, 1986, Whittle 1985b; 
Appenzeller & Ostreicher 1988; De Robertis & Shaw 1990; 
Veilleux 1991b). 

By comparing the FWHM of [O i] 26300 to that of [S n] 
>16716 in Seyfert galaxies, and by comparing the FWHM of 
[O hi] 24363 to that of [O hi] 25007 in Seyfert 2 galaxies (in 
which there is little or no contamination from broad Hy 
and/or Fe n emission), we have argued that the correlations 
are best interpreted as fundamentally with ionization poten- 
tial rather than with critical density (Moore & Cohen 1994, 
hereafter MC). We also proposed a new interpretation that 
these correlations are a column density effect, in which case 
the correlations are also consistent both with constant ion- 
ization parameter and with velocity increasing outward. As 
we pointed out in MC, line width is correlated with neither 
critical density nor ionization potential in some objects. In 
these objects, the similarity between profiles spanning a 
broad range of ionization potential and critical density is 
inconsistent with a scenario in which lines of low and high 
IP/w cr are formed in physically distinct regions. The [Fe vii] 
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26087 and [Fe x] 26374 lines, which are very sensitive to 
ionization parameter (MC), are particularly important in 
this regard, and we emphasize these lines here. 

In § 3, we present a sample of narrow-line profiles in 
intermediate Seyfert galaxies. The sample was chosen with 
the goal of obtaining high-quality [Fe vii] 26087 and, when 
possible, [Fe x] 26374 profiles. The sample consists of 
MCG 8-1 1-1 L Mrk 79, Mrk 704, Mrk 841, NGC 4151, and 
NGC 5548. 

We discuss briefly details of the observations and of the 
data reduction and profile extraction in § 2. In § 3, we 
discuss general implications for the NLR of line profile fea- 
tures. This discussion is based on previous observational 
studies as well as on our own. It will be useful to have 
already highlighted those general features of NLR spectra 
that discriminate best between different NLR models when 
we discuss our spectra and profiles. In § 4 we discuss specific 
implications of our data for NLR models. This discussion is 
facilitated by references to models presented in a compan- 
ion paper (Moore & Cohen 1996, hereafter Paper II). The 
most constraining feature of the observed line profiles is the 
similarity in some objects of the profiles of lines correspond- 
ing to a wide range of IP/fl cr , and with very different sensiti- 
vities to ionization parameter. In § 5, we confirm for our 
sample the well-known correlations of velocity shift and 
width with IP and n CT . However, our sample, which includes 
one contradictory object (NGC 5548), is too small to 
confirm by itself the conclusion of MC that the correlations 
are fundamentally with ionization potential rather than 
with critical density. Finally, we summarize our results in 
§ 6 . 

2. DATA ANALYSIS 
2. 1 . Observations and Reductions 

Observations were made with the 200 inch (5.08 m) Hale 
telescope and double spectrograph (Oke & Gunn 1982) on 
six nights in 1983 and on a seventh night in 1984. The 
double spectrograph incorporates a dichroic beam splitter 
to cover approximately 3100-5200 A on the blue side and 
5200-11000 A on the red side. The blue side had an RCA 
320 x 512 CCD with 30 fim pixels, which corresponds to a 
spatial scale of 078 pixel - 1 . One DN (count) corresponds to 
6 photoelectrons, and the readout noise was 8 DN pixel” 1 . 
A 1200 g mm' 1 grating blazed at 4700 A yields a dispersion 
of 1.09 A pixel - 1 at 3500 A (200 km s" 1 FWHM). The red 
side had a TI 800 x 800 CD with 15 fim pixels, which corre- 
sponds to a spatial scale of 076 pixel -1 . One DN corre- 
sponds to 2 photoelectrons, and the readout noise was 5 
DN pixel 1 . A 1 200 g mm " 1 grating blazed at 7 1 00 A yields 
a dispersion of 0.82 A pixel 1 at 6300 A (75 km s _1 
FWHM). 

Because of the small wavelength coverage, most observ- 
ations required unique wavelength settings. Grating tilts 
were set by hand with a vernier reading, which could not be 
reset to a specified wavelength with high precision. For 
most observations, wavelength calibration lamps and flat- 
field exposures were taken immediately before or after the 
observation at the hour angle and declination of the start or 
end of the observation, in order to minimize errors due to 
flexure in the spectrograph. Occasionally, calibration expo- 
sures were taken between two identical setups on an object. 
All observations were made with a 1" slit except one (24363 
in NGC 5548). 


For all the data reduction, we used the astronomical data 
reduction software package VISTA, developed at Lick 
Observatory. Straightforward procedures were used for 
baseline and LED preflash subtraction (for exposures in the 
blue) as well as for flat-fielding. 

For each column, the spatial profile of the sky back- 
ground was modeled by a fourth-order polynomial fit 
outside those rows bracketing the spatial profile of the 
object. The procedure includes automatic detection and 
rejection of cosmic-ray events as well as of hot and cold 
pixels in the sky background. From each column of the 
exposure, we subtracted the sky background model for that 
column. 

The orientation of the CCD was skewed slightly with 
respect to the grating rulings, resulting in a shift in the 
centroid of the spatial profile of about 3 pixels over the 
length of the spectrum. The spectrum was extracted from a 
symmetric window, typically 6 8 pixels wide, centered on a 
parabolic fit to the centroid curve. The extraction included 
interpolation to allow for fractional pixels. 

The size of the window used for spectral extractions was 
chosen to be roughly the apparent full width at zero- 
intensity of the spatial profile of the nuclear continuum. 
Observations were acquired under a range of atmospheric 
seeing conditions as follows. Spectra of MCG 8-11-11 were 
acquired under 3" -4" seeing. Spectra of Mrk 79 were 
acquired under 1'.'5 seeing. Separate spectra of Mrk 704 
were acquired under 1", F.'5, 2", and 4"- 5" seeing. Spectra of 
Mrk 841 were acquired under 175 and 3" seeing. Separate 
spectra, covering the same wavelength region, of NGC 4151 
were acquired under F.'5 (Fig. 5a), and 4"-5"(Fig. 5b) seeing. 
Spectra of NGC 5548 were acquired under 175-2" seeing. 
Even for NGC 4151, the nearest Seyfert galaxy in the 
sample, and for which separate spectra were acquired under 
the most extreme differences in seeing, spectra obtained 
under different seeing conditions showed very small, if any, 
systematic differences. 

The extracted spectra were wavelength calibrated using 
polynomial fits to the measured centroids and known wave- 
lengths of approximately 10 comparison lines in the blue 
and approximately 15 lines in the red. The order of the 
polynomial ranged between 2 and 4 and resulted in an rms 
residual of 0.02-0.06 A. The absolute accuracy of the wave- 
length calibration can be inferred only from the night-sky 
lines in the spectrum. Offsets between the wavelength- 
calibrated object spectra and those of the night sky were of 
the order of 0.2 A in scans with night-sky lines strong 
enough to permit a reliable comparison. The absolute accu- 
racy could be better than 0.2 A because this is roughly the 
precision with which the centroids of night-sky lines can be 
calculated. After extracting each spectrum and correcting 
for atmospheric extinction, the spectrum of a standard star 
was used to flux calibrate the object spectrum. 

Occasionally an emission line of interest is close in wave- 
length to an atmospheric absorption feature. To correct for 
the feature, we used a transmission curve derived from a 
normalized standard star spectrum, in which the depth of 
the absorption was determined empirically by comparing 
the depths of the absorption troughs in the object and stan- 
dard star (which were observed at slightly different air 
mass). 

There were two profiles that were seriously corrupted by 
atmospheric absorption, those of [S n] 26731 in MCG 8-11- 
1 1 and of [S n] 26716 in Mrk 79 (both by B band). We have 
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excluded these profiles from formal presentation and 
analysis. However, profile parameters were calculated none- 
theless, and after correcting for the absorption, the param- 
eters calculated for 26716 in MCG 8-11-11 were in 
reasonable agreement with those calculated for 26731. Simi- 
larly, the parameters of 26731 in Mrk 79 agreed with those 
of 26716. In Mrk 704, the B-band feature ends just blueward 
of [S n] 26716, and inspection of the standard star shows 
that the atmospheric absorption does not affect the emis- 
sion line directly. Nevertheless, the absorption complicated 
selection of the blue wing of the 26716 profile enough that, 
after subtracting off an estimated underlying red wing of Ha 
broad, we considered the profile unsuitable for further 
analysis. The atmospheric absorption feature at 6280 A is 
just redward of the [Fe vii] 26087 profile in Mrk 704, but 
again inspection of the standard star spectrum shows that 
the line profile is not alfected directly. The [Fe vii] line of 
Mrk 704 was not excluded from further analysis because the 
strength and shape of the adjacent continuum permitted 
reliable subtraction, despite the absorption feature. The [Ne 
hi] 23869 profile in MCG 8-11-11 is just blueward of Ca n 
absorption at 3935 A from the Galaxy, but the absorption 
feature is narrow, and again neither the profile nor the con- 
tinuum subtraction are affected seriously. 

Several of the emission lines were observed in more than 
one scan either because several observations were made 
with the same grating tilt or because spectra overlapped 
enough in wavelength to include a common emission line. 
Such scans were averaged weighting by (S/N) 2 , where the 
signal-to-noise ratio (S/N) in each spectrum was determined 
empirically from the rms scatter about a low-order poly- 
nomial fit to a featureless part of the continuum. 

2.2. Profile Extraction 

We chose the limits of each emission line by visual inspec- 
tion and subtracted a polynomial fit to a featureless part of 
the adjacent continuum. The profiles were then placed on a 
velocity scale. Small heliocentric corrections were applied to 
the velocities as appropriate. In the ideal case, the zero 
point of the velocity scale corresponds to the recessional 
velocity of the NLR center of mass. The least subjective 
estimate of the NLR center-of-mass velocity is that of 21 cm 
radiation. We used a 21 cm measurement, taken from 
Heckman, Balick, & Sullivan (1978), for Mrk 79, NGC 
4151, and NGC 5548. A 21 cm measurement was not avail- 
able for the remaining objects. We used the peak of [O in] 
25007 for Mrk 704, the centroid of [O i] 26300 for MCG 
8-11-11, and the peak of [O i] 26300 for Mrk 841. The 
appropriate measurement was picked from observational 
considerations such as the S/N, profile asymmetry, and/or 
the degree to which a profile peak could be determined 
unambiguously, rather than from model-dependent con- 
siderations in which a particular measurement is assumed a 
priori to best represent the galactic systemic velocity. 

In some cases, the spread in cloud velocity is large 
enough to blend together the individual profiles of two 
nearby emission lines. To reconstruct the individual pro- 
files, we used a deblending program developed at UCSD by 
Steve Bradley. The program allows a template profile to be 
defined from the data, from which model profiles can be 
derived by rescaling the template in intensity and by using 
interpolated profile fluxes to vary the wavelength and width 
of the template. 

Further details of the appropriate deblending procedure 


depend on the nature of the blend. We describe the deblend- 
ing of four particularly important lines, [O m] 24363 (high- 
density diagnostic) from broad and narrow Hy, [Fe x] 
26374 (high-ionization parameter diagnostic) from [O i] 
26363, [O i] 26300 (large column density diagnostic) from 
[S hi] 26312, and [S ii] 226716, 6731 (low-density 
diagnostic). 

To subtract Hy from [O in] 24363, we used the broad 
and narrow components of Up to model the respective 
components of Hy, scaling them in intensity and width and 
varying the position until the final profile matched, in a 
least-squares sense, that portion of the Hy profile not con- 
taminated by [O in] 24363 emission. We did not assume 
that the ratios of the narrow and broad components are the 
same for both Hy and H/?. For Mrk 79 and NGC 5548, an 
analogous procedure was used to deblend [S in] 26312 
from [O i] 26300 using [O hi] 25007 as a template [O i] 
profile. The net result here is to replace that portion of the 
[O i] profile contaminated by [S in] emission by a template 
[O in] 25007 profile, where the template is fitted in intensity 
and width to match that portion of the [O i] profile not 
contaminated by [S m] emission. (The [O i] and [O hi] 
profiles are similar to each other in Mrk 79 and NGC 5548.) 
To deblend [S hi] from [O i] in NGCC 4151, we shifted and 
scaled the [O i] + [S in] blend itself to model the [S m] 
profile and subtracted the result from the [O i] + [S in] 
blend. To subtract [O i] 26363 from [Fe x] 26374, we used 
the profile of [O i] 26300 scaled by 0.33, where the scale 
factor is derived from atomic parameters. To deblend [S n] 
226716, 6731, we used an iterative procedure with the 
blended 26731 profile as an initial template. We fitted the 
blended 26731 profile to that of the blended 26716 profile 
with an intensity scale factor as the single free parameter. 
The resulting profile is a model of a deblended 26716 
profile. The model is poorly determined because the tem- 
plate 26731 profile is highly contaminated by 26716 emis- 
sion. Nevertheless, the model 26716 profile is subtracted 
from the blend at 26716, resulting in a partially deblended 
26731 profile. The procedure is then repeated using the par- 
tially deblended 26731 profile as a template. The result is a 
better model of a deblended 26716 profile than that of the 
previous iteration because the 26731 template is less con- 
taminated by 26716 emission. More than three iterations 
only degraded the S/N with little or no improvement in the 
deblending. 

We did not apply a formal deblending procedure to 
reconstruct profiles of [N n] 26548, Ha (narrow 
component), and [N n] 26584. The only objects for which 
we attempted to analyze these profiles were those in which 
the narrow components were not blended with each other 
and the underlying component of broad Ha could be dis- 
cerned by visual inspection and subtracted off in much the 
same way we subtracted off the continuum in other lines. 
The H/? profile was used merely as a guide for the 
44 continuum placement.” 

An atlas of the line profiles extracted from the reduced 
spectra is given in Figures 1-6. All the profiles presented 
here are normalized to the peak. Lower resolution studies, 
in which a large spectral range was obtained on a single 
night, resulted in more reliable relative flux estimates than 
did our own study. We include these relative flux estimates 
in Table 1. The dereddened emission-line fluxes for MCG 
8-11-11, Mrk 79, Mrk 704, Mrk 841, and NGC 5548 are 
from Cohen (1983), and those (not corrected for reddening) 
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TABLE 1 


Narrow Emission-Line Flux Relative to H0 


Line 

MCG 8-11-11 

Mrk 79 

Mrk 704 

Mrk 841 

NGC 4151 

NGC 5548 

[O ii] A3727 

3.24 

2.20 

0.66 

0.99 

1.47 

1.33 

[Fe vn] 23760 

0.21 

0.25 

1.06 

0.47 

0.14 

0.52 

[Ne hi] 23869 

1.46 

1.57 

1.25 

1.62 

1.35 

1.88 

[O m] 24363 

0.32 

0.48 

0.53 

0.18 

0.34 

1.23 

He ii 24686 

0.16 

0.25 

0.43 

0.32 

0.15 

0.26 

H i 24861 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

[O hi] 25007 

8.46 

10.11 

6.45 

11.08 

15.29 

9.61 

[Fe vn] 25721 

0.08 

0.17 

0.52 

0.07 

0.11 

0.10 

He i 25876 

0.10 

0.13 

0.11 

0.17 

0.15 

0.07 

[Fe vii] 26087 

0.10 

0.28 

0.74 

0.30 

0.22 

0.35 

[O i] 26300 

0.40 

0.36 

0.16 

0.32 

1.06 

0.26 

[S hi] 26312 


0.05 

0.12 

0.04 

0.03 

0.03 

[Fe x] 26374 

0.06 

0.11 

0.49 

0.11 

0.12 

0.08 

H i 26563 

2.08 

2.93 

2.62 

2.87 

5.29 

2.43 

[N ii] 26584 

2.06 

2.16 

0.91 

1.08 

2.24 

1.10 

[S ii] 26716 

0.58 

0.72 

0.22 

0.35 

0.82 

0.21 

[S ii] 26731 


0.72 

0.31 

0.37 

1.00 

0.20 


for NGC 4151 are from Boksenberg et al. (1975). The [Fe x] 
26374 line strength in NGC 4151 was determined from the 
26374 to 26300 line strength ratio of Veilleux (1988) and the 
26300 line strength of Boksenberg et al. (1975), as 26374 was 
not included in the latter study. These line strengths have 
implications for some of the models discussed in Paper II. 
For possible future reference, we have also included [S in] 
26312 relative line strengths based on our measured [S hi] 
26312 to [O i] 26300 line strength ratios and the previously 
published [O i] line strengths. 

3. GENERAL IMPLICATIONS FOR THE NLR 

Before discussing interpretations of the line profiles, we 
summarize some important characteristics of NLR model 
profiles that should not be overlooked. The basic goal is to 
provide interpretations that stem directly from the observed 
profiles to the extent possible, with only minimal reference 
to specific models. The most important premise is that the 
preponderance of blueward asymmetry constrains the 
highest velocity clouds to move mostly in the radial direc- 
tion and precludes the presence of a broad distribution of 
velocities at each radial coordinate, r (asymmetric distribu- 
tions of gas would produce red and blue asymmetries with 
equal likelihood). Neglecting obscuration, a spherical shell 
of delta function-emitting clouds moving radially with 
velocity v 0 yields a boxcar profile with half-width t? 0 . As can 
be inferred directly by visual inspection, replicating the 
observed profiles requires that either boxcar profiles with a 
range of velocity width must be superposed, or else assump- 
tions such as purely radial cloud flow or a spherical NLR 
must be relaxed. 

In MC we showed that the profile wings are not sensitive 
to either the orientation or the opening angle of a biconical 
line-emitting structure (with the caveat that the velocity 
field is renormalized appropriately when the line of sight 
falls outside the biconical structure), even though the profile 
core is. In Moore (1994, 1995), it is demonstrated that the 
blue wing is not very sensitive to details of the obscuration, 
whether internal or external to the emitting clouds. Thus, it 
is possible to interpret the blue wings readily in the context 
of a simplified model NLR even though collimated emission 
and obscuration are likely present. This point is worth 
emphasizing because it follows then from the absence of a 
broad distribution of cloud velocities at each r (see above) 


that the range in height and in velocity width of boxcar 
profiles that must be superposed to replicate the observed 
blue wings relates directly to the spatial dependence of the 
line emission. For lines that are useful diagnostics of elec- 
tron density (n e \ column density (jV co1 ), and/or ionization 
parameter (T), the spatial dependence of the line emission in 
turn corresponds directly to how n e , N col , and T vary spa- 
tially throughout the NLR. The effects of spatial variations 
in n e , N coh and T are discussed qualitatively when we 
discuss the profiles of individual objects. The effects are 
illustrated more explicitly (though less independent of 
detailed modeling) in Paper II. 

It is an observational fact that, in some objects, line pro- 
files of high and low IP have very similar blue wings (see the 
discussion in MC as well as the profiles presented in § 3 and 
in Veilleux 1991a). It follows that if the velocity field is at all 
close to monotonic, then there is emission into lines of both 
high and low IP at each /*, which could arise in either of two 
scenarios. 

The first possible scenario is that individual clouds are 
stratified in ionization. The second possible scenario is that 
there is a broad distribution in the electron density of 
unstratified clouds and, therefore, in T at each r. Through- 
out this paper, we use the term “stratified cloud " to denote 
a single cloud of sufficient ionization stratification to emit 
[Fe x] and/or [Fe vn] lines from the front as well as [O i] 
and [S ii] lines from the back. An “unstratified cloud" is 
optically thin to the photoionizing continuum so that a 
single ionization state describes the entire cloud. It is worth 
pointing out that blueward asymmetry and stratified clouds 
are consistent with either inflow or outflow according to 
whether the obscuration is mostly internal or external, 
respectively. Because the applicability of the first scenario is 
fundamental to our detailed modeling (Paper II) and to 
several of our conclusions, we outline here why we believe 
the first scenario applies. We note, however, that work 
favoring the second scenario is underway (Komossa & 
Schulz 1994; Tsvetanov et al. 1994). Our arguments for 
stratified clouds are as follows : 

1. Clouds of very different n e , but at the same r, cannot 
all be pressure confined. Because the equilibrium tem- 
perature of photoionized gas varies only slightly with n e 
(MC), there is a narrow range of n e at each r that can be in 
pressure equilibrium either with the intercloud medium 
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(ICM) or with a magnetic field. Because the cloud expan- 
sion timescale is a small fraction of the crossing time, a 
broad distribution of n e at each r would require that clouds 
are continually created and destroyed (this has been sug- 
gested for broad-line region [BLR] clouds based on a sta- 
bility analysis [Krinsky & Puetter 1992]). 

2. A contrived distribution of n e must be invoked for the 
second scenario so that each n e contributes comparable flux 
to the appropriate emission line. Considering that the 
volume emissivity for optically thin gas is j oc n], and the 
emission per fixed mass is proportional to n e , arranging for 
clouds differing by several decades in n e to produce compa- 
rable emission-line flux at each r, as observations require, is 
not trivial. 

3. An inflow model requires stratified clouds to explain 
the preponderance of blueward asymmetry. An outflow 
model (with external obscuration) can accommodate 
unstratified clouds, but a mechanism must be found to out- 
wardly accelerate clouds of very different density to nearly 
the same velocity. Otherwise, a broad distribution of pro- 
jected velocity at each r, and consequently symmetric pro- 
files, would result. In principal, clouds of very different n e 
can be accelerated to similar velocities over very different 
length scales. This scenario could produce distinct 
coexisting velocity fields, with different lines coming from 
different underlying cloud populations at different r. 
However, to account for the striking similarity in the pro- 
files of different lines in some objects, this exceedingly con- 
trived scenario requires that the variation of emission-line 
flux with projected velocity is virtually identical for lines in 
which the variation of flux with r is very different. 

4. While lines of high and low IP/n cr are very similar in 
some objects, the high IP/n cr lines are systematically 
broader in other objects. Both types of objects would not be 
explained naturally by a broad distribution in n e and T at 
each r. Actually, both types of objects are also not easily 
explained in the classical picture of the NLR in which the 
highest velocity clouds are at small r , where lines of low n cr 
are collisionally suppressed, resulting in systematically nar- 
rower low n cr lines. To explain both types of objects in the 
classical picture of the NLR, the two types of objects must 
differ in the maximum n e by several decades, which is the 
relevant range of n cr . In Paper II, we present models in 
which line width correlations are a natural consequence of 
outward acceleration and of N co , oc r -4/3 , the latter 
resulting from spherical clouds of constant mass and with 
constant ionization parameter. In these models, objects in 
which the low-ionization potential lines are broad and 
narrow differ in AT col by factors of only order unity. 

4. INDIVIDUAL OBJECTS 
4.1. MCG 8-11-11 

The profiles extracted from MCG 8-11-11 spectra are 
plotted in Figure 1. We used the [O i] 26300 line rather 
than [O hi] 25007 to estimate the systemic velocity of this 
object because 26300 is more symmetric than 25007, and an 
asymmetric profile yields an underestimated systemic veloc- 
ity due to obscuration effects. We used the centroid rather 
than the peak because the peak of [O i] 26300 is fairly 
flat-topped, so the peak location is ambiguous. 

The [O n] 23727 profile is a composite of the doublet, 
23726 + 23729, corresponding to a separation of 217 km 
s“ *. The blue wing of 23727 arises mostly from 23726 emis- 


sion. Clearly, from the apparent extent of the blue wing in 
Figure la, there are clouds with v as — 1000 kms' 1 appre- 
ciably emitting 23726. The critical density of this line is 
3 x 10 3 cm -3 , so the electron density of the v « — 1000 km 
s -1 clouds emitting 23726 does not greatly exceed 10 4 * 
cm -3 . 

As is the case with the profiles of Veilleux (1991a), the 
presence of a substantial velocity gradient across the NLR 
can be inferred by visual inspection of the profiles. In partic- 
ular, note that a similar range of radial velocity can be 
inferred for all the profiles shown in Figure 1 (i.e., for each of 
the profiles, a substantial range in velocity width of boxcar 
profiles must be superposed to replicate the blue wing). In 
the absence of a broad distribution of cloud velocity at each 
r (§ 3), the similar range of radial velocity corresponds to 
similar spatial extents for the line emission. In the case of 
[Fe vn] 26087, the similar spatial extent suggests a roughly 
constant ionization parameter throughout the NLR 
because [Fe vii] 26087 is much more sensitive to ionization 
parameter than the other lines (MC). In the case of [O i] 
26300, the similar spatial extent suggests that the clouds are 
also optically thick to the photoionizing continuum 
throughout the NLR. The same consideration could also 
suggest that the low-density limit ( n e n cr ) is valid through- 

out the NLR for even the lowest critical density lines ([O n] 
23727 and [S ii] 226716, 6731) because the effects of col- 
lisional de-excitation are apparent in neither the wings nor 
in the cores of these lines. Although the effects of collisional 
de-excitation on the profiles are illustrated explicitly in the 
detailed kinematical models of Paper II, they can also be 
discerned qualitatively as follows. If the NLR is spherical, 
and if obscuration can be neglected, and if the velocity field 
is radial even for the lowest velocity clouds, then the line 
profile is a simple superposition of boxcar profiles. Col- 
lisional de-excitation effects in the highest velocity clouds 
would be evident if the low critical density line profiles were 
lacking a contribution from the highest velocity boxcar pro- 
files (i.e., the profile wings). Collisional de-excitation effects 
in the lowest velocity clouds would be evidenced by a broad 
flat profile peak (there would be no contribution from the 
lowest velocity boxcar profiles). Absence of collisional de- 
excitation effects throughout the [O n] 23727 and [S ii] 
226716, 6731 profiles would then limit the cloud electron 
density to n e <, 10 4 cm -3 throughout the NLR. However, 
the NLR model we favor includes n e « 10 6 cm -3 gas (see 
Paper II, MC, and Moore 1994) and requires emission from 
low-density and low-velocity gas outside the (parameter 
space of the) NLR model to account for the [O n] 23727 
and [S ii] 226716, 6731 profile cores. 

4.2. Mrk 19 

The profiles extracted from Mrk 79 spectra are plotted in 
Figure 2. We have presented the profile of [S ii] 26731 
rather than that of [S ii] 26716 because the latter is badly 
corrupted by the B-band absorption feature. The [Fe x] 
26374 profile is of such poor quality that almost nothing 
can be said about the NLR based on this profile except that 
there are [Fe x] 26374-emitting clouds with v « — 1000 km 
s’ 1 and that there is some range of (radial) velocity sampled 
by the 26374-emitting clouds because the profile is not a 
single boxcar profile. However, 26374 is so sensitive to ion- 
ization parameter (MC) that even this small piece of infor- 
mation is useful and implies that the ionization parameter is 
roughly constant over whatever range of velocity can be 
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Fig. 1. — Narrow emission-line profiles extracted from spectra of MCG 8-11-11. A scenario in which the NLR can be decomposed into physically distinct 
regions in which the low-ionization potential/critical density lines and the high-ionization potential/critical density lines are formed is inconsistent with the 
similarity in the profiles shown here. To account for the similarity in the profiles requires significant emission in each line over a large fraction of the NLR. In 
particular, note the similar velocity ranges sampled by [S n] A6716 and [O i] ,16300. These lines differ by three decades in critical density. Panels {g) and (h) 
represent [S n] before and after deblending, respectively. The trough of the B-band absorption feature coincides with the peak of A6731. Due to the 
coincidence in wavelength, the B-band correction is poor, as evidenced by the residual feature at the peak of A6731. Although A6731 was excluded from 
further analysis, profile parameters were calculated, nonetheless, and they were in reasonable agreement with those of A6716. 


inferred from the profile. Note that the [Fe vn] A6087 
profile has an extended blue wing that is not present in 
other profiles (we exclude the [Fe x] A6374 profile from the 
discussion owing to the very low S/N in the wings). The 
FWHM of [Fe vn] A6087 is similar to that of [O in] /15007, 
but a width parameter defined at a lower fraction of the 
peak intensity is significantly larger for [Fe vii] than for 


[O hi]. Our favored models (n e ocr~ 2 , veer 112 , and 
N col oc r“ 4 ' 3 ) show the same effect as well(Paper II). 

In both objects, MCG 8-11-11 and Mrk 79, the [Fe vn] 
A6087 profile is broad. The [O i] A6300 profile is either 
broad (MCG 8-11-11) or narrow (Mrk 79). When [O i] is 
broad, [S n] is also broad, and when [O i] is narrow, [S ii] 
is also narrow. This is consistent with the trend discussed 
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Fig. 2.— Narrow emission-line profiles extracted from spectra of Mrk 79. Note that in Mrk 79, unlike in MCG 8-11-11, [O i] and [S n] are both narow 
relative to [Fe vii]. We explain the difference between MCG 8-11-11 and Mrk 79 in terms of identical NLRs except that clouds at the inner radius of MCG 
8-11-11 start out with somewhat higher column density than those of Mrk 79. In this scenario for Mrk 79, 7V col /7V inifl decreases to less than one before the 
clouds reach their maximum velocity at the outer radius, while in the MCG 8-11-11 scenario, N <JOl /N min > 1 throughout the NLR. 


in MC for a much larger sample (^50) of objects. The 
trend for [O i] and [S ii] to have similar widths poses a 
difficulty for interpreting line width correlations as funda- 
mentally with n cr because these lines differ in n cr by three 
decades. 

We can explain the difference between MCG 8-11-1 1 and 
Mrk 79 in terms of identical NLRs, except that clouds at the 
innermost radius in MCG 8-11-11 start out with somewhat 
higher column density than those in Mrk 79. In Mrk 79, 
N col decreases to less than N min before the clouds reach their 


maximum velocity at the outermost radius, while in MCG 
8-11-11, N col /N min > 1 throughout the NLR. We define 
N min as the minimum column density needed for an [O i] 
26300-emitting zone to exist in the back of the clouds. 
Referring to the models of Paper II, MCG 8-11-11 corre- 
sponds to model 3 and Mrk 79 corresponds to model 4. We 
postulate that our model for Mrk 79 is representative of 
objects that show correlations of line width with ionization 
potential, while our model for MCG 8-11-11 is representa- 
tive of objects that do not. 
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4.3. Mrk704 

The profiles extracted from Mrk 704 spectra are plotted 
in Figure 3. The systemic velocity was estimated from the 
peak of [O hi] A5007. Note the extended blue wing in the 
[Fe vii] profile, from which it is clear that the [Fe vn]- 
emitting clouds sample a wide range of velocity. Of even 
greater significance, a substantial range of velocity can also 
be inferred from the [Fe x] ,16374 profile, which is extremely 
sensitive to ionization parameter (MC). Because of the wide 
range of cloud velocity that can be inferred directly from the 
observed [Fe x] A6374 profile, our case for constant ioniza- 



- 1 000 0 1 000 


—I 1 1 * 1 1 r 

(e) MRK 704 : [Fe Vll] \6087 



I 1 1 1 ■ 

-1000 0 1000 


(g) MRK 704 : [Fe X] A6374 



-1000 0 1000 
velocity (km/sec) 


tion parameter is strongest for Mrk 704, NGC 4151, and 
NGC 5548. 

As is the case for Mrk 79, the extended blue wing in the 
[Fe vii] profile of Mrk 704 is not evident in the profiles of 
[O i] >16300 and [S n] >U6716, 6731. The obviously low 
quality of the [O i] and [S n] profiles precludes a more 
definitive statement, however. The [O i] and [S ii] profiles 
of Mrk 79 are suggestive of the following interesting effect, 
which is evident in model 4 of Paper II though the effect is 
inconclusive here owing to the low quality of these profiles. 
The effect is that a width parameter defined in the wings is 


—i 1 » 1 r i 

(b) MRK 704 : [Fe Vll] A3760 



-1000 0 1000 



-1000 0 1000 



velocity (km/sec) 


Fig. 3. — Narrow emission-line profiles extracted from spectra of Mrk 704. As is the case in Mrk 79, the extended blue wing of the [Fe vii] profile of Mrk 
704 is not evident in the [O l] and [S n] profiles 
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larger for [Fe vii] than for [O i] and [S n], even though the 
FWHM is about the same, and the width at any fraction of 
the peak intensity is similar for [O i] 26300 and [S ii] 26716. 

4.4. Mrk 841 

The profiles extracted from Mrk 841 spectra are plotted 
in Figure 4. The systemic velocity was estimated from the 
peak of [O i] 26300. Note the very extended (>1000 km 
s _1 ) blue wing of the [Fe vii] 26087 profile. This wing is not 
present in other profiles (the [Fe x] profile is inconclusive). 


For Mrk 841, unlike Mrk 79 and Mrk 704, the width at any 
fraction of the peak intensity is greater, almost up to the 
peak itself, for [Fe vii] than for [O i] and [S n]. Note also 
that the [S n] profiles are narrow, and the [O i] profile is 
narrow as well, in keeping with the trend discussed here and 
in MC. Thus, we assert that Mrk 841 is best explained by 
model 4. 

4.5. NGC4151 

The profiles extracted from NGC 4151 spectra are 
plotted in Figure 5. We include two profiles (Figs. 5a~5b) of 





Fig. 4. — Narrow emission-line profiles extracted from spectra of Mrk 841. Note that [O i] and [S n] are both very narrow relative to [Fe vii]. This is in 
keeping with a general trend: when [O i] is narrow, [S ii] is narrow as well, and when [O i] is broad, [S n] is broad also. These lines differ in critical density 
by three decades, however. 
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Fig. 5. — Narrow emission-line profiles extracted from spectra of NCG 4151. We include two profiles, (a) and (b), of [O i], [Fe vn], and [Fe x] 
corresponding to two observations 105 days apart because the Fe lines of NGC 4151 have been reported in the literature as possibly variable (Penston et al. 
1984). Note that the [Fe vii] profiles of NGC 4151 and Mrk 704 are virtually identical, suggesting a very simple common acceleration mechanism. Note also 
that the [Fe x] profile has the same overall characteristics as the [Fe vn] profile of any object presented thus far. This would not likely result if [Fe x] and 
[Fe vn] were emitted in physically distinct regions. 


[O i] A6300, [Fe vn] ,16087, and [Fe x] A6374 correspond- 
ing to two observations 105 days apart. NGC 4151 is highly 
variable in its continuum luminosity, and the forbidden Fe 
lines of NGC 4151 have been reported as possibly variable 
(Penston et al. 1984). We see no compelling evidence of 
variability in the profiles here (of course, our time baseline is 
much smaller than the 3 year baseline of Penston et al.). 

As in every object presented thus far, the [Fe vn] profile 
of NGC 4151 shows a very extended (> 1000 km s -1 ) blue 
wing. The similarity in the [Fe vn] profiles of NGC 4151 


(Fig. 5a) and of Mrk 704 is striking. The profiles are vir- 
tually identical. In all five objects whose profiles we have 
presented thus far, the [Fe vn] profiles are very similar and 
do not differ even in an overall scale factor in the velocity 
width. In every object, we can infer [Fe vn] emission out to 
v % 1000 km s" 1 and not much higher. The remarkable 
similarity in the [Fe vii] profiles of each object suggests a 
very simple common acceleration mechanism. Likewise, 
Whittle (1991) has pointed out previously that the similarity 
in the [O hi] A5007 profiles from a large number of objects 
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Fig. 6. — Narrow emission-line profiles extracted from spectra of NGC 5548. Note that the [O i] profile is broad, but the [S n] profile is narrow. This is in 
contrast to a general trend discussed here and in MC. 

suggests a common acceleration mechanism for the 25007- 4151 are more representative of a typical Seyfert 1 galaxy 

emitting gas. than those of Mrk 704. We were able to extract reliable 

The similarity in the [Fe vn] 26087 profiles of NGC 4151 26087 and 26374 profiles of NGC 4151 because our spectra 

and Mrk 704 has possible implications concerning how rep- have high S/N due to this object’s brightness. Of the six 

resentative our sample of Seyfert 1 galaxies is. Because the objects in our sample, the 26087 line of Mrk 704 is the 

sample was chosen with the goal of obtaining high-quality strongest relative to Hj? (0.87), while that of NGC 4151 is 
[Fe vii] 26087 and, when possible, [Fe x] 26374 profiles, the weakest (0.07). The object in our sample with the strong- 

selection was biased toward objects with strong 26087 and est 26087 line and an object that is representative in its 

26374 emission. It is quite possible that the same properties 26087 line strength of a typical Seyfert 1 galaxy have vir- 

of the same that allow us to obtain reliable 26087 and 26374 tually identical profiles, which suggests that the density and 
profiles also render the sample unrepresentative of Seyfert 1 velocity profiles throughout the NLR are very similar in 

galaxies in general. The forbidden Fe line strengths of NGC both objects. That is, we do not infer any important differ- 
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ences between Mrk 704 and NGC 4151 other than the dif- 
ferent [Fe vii] line strengths. The argument could be tested 
by selecting at random one or two objects with weak 
[Fe vii] line strengths and investing the telescope time 
required to obtain S/N comparable to that of our profiles of 
[Fe vii] 26087 in both Mrk 704 and NGC 4151. We note 
also that, although Veilleux (1991a) did not include any 
[Fe x] 26374 profiles and included [Fe vii] 26087 profiles 
for only Mrk 359 and NGC 4151, the lower ionization line 
profiles in his sample are similar to those in our own 
sample. 

Any of the [Fe vii] profiles presented thus far will suffice 
more or less as an archetypical [Fe vii] 26087 profile. The 
[Fe x] 26374 profile of NGC 4151 has the same overall 
characteristics as the [Fe vii] 26087 profile of any object. 
This suggests that [Fe x] is emitted in a region that is 
similar spatially to the region of [Fe vii] emission rather 
than from a physically distinct region close to the central 
source with n e « n cr = 5 x 10 9 cm -3 . Our case for 
constant-ionization parameter, based on the forbidden Fe 
lines, is probably strongest for NGC 4151. In Paper II we 
present some models in which ionization parameter 
decreases outward but not so much that the model [Fe vii] 
26087 profile is inconsistent with observations. The [Fe x] 
26374 profile, on the other hand, is extremely sensitive to 
ionization parameter. Almost any spatial variation at all in 
the ionization parameter results in a synthetic [Fe x] 26374 
profile that is, contrary to observations, different from that 
of [Fe vii] 26087 and other profiles. 

NGC 4151 is like MCG 8-11-11 in that the [Fe vii] and 
[O i] profiles are very similar. Thus, we infer similar spatial 
extents to the [O i] and [Fe vii] line emission, which 
requires that either [O i] and [Fe vii] are emitted from 
clouds at the same r but with very different n e , or else the 
emitting clouds are individually stratified in ionization. We 
argued for the latter in § 3. We do not have a spectrum of 
NGC 4151 that includes the [S n] lines, but the blue wing of 
the [S ii] 26716 profile presented by Veilleux (1991a) is very 
similar to that of [O i] 26300 (either ours or Veilleux’s). 
Again, this agrees with the trend discussed here and in MC, 
namely, when [O i] is broad, [S ii] is broad as well. 

4.6. NGC 5548 

The profiles extracted from NGC 5548 spectra are 
plotted in Figure 6. The overall characteristics of the 
[Fe vii] profile from NGC 5548 are consistent with the 
[Fe vii] profile from the above five objects, although we 
cannot infer as wide a range of cloud velocity from this 
profile. The latter could be due to a smaller v ma Jv min for the 
[Fe vii] 26087-emitting gas in this object or perhaps to a 
larger chaotic component of the velocity field. We define 
i; max as the maximum radial velocity of NLR clouds. It is the 
velocity at the innermost or outermost NLR radius, accord- 
ing to whether the cloud velocity increases inward or 
outward. Similarly, v min is the minimum radial velocity of 
NLR clouds. When the ratio of the maximum radial veloc- 
ity in the NLR to the minimum radial velocity (tWxMmn) is 
decreased, the wings are less pronounced, and the result is a 
profile closer to that of a boxcar function. 

Note that our discussion of the [Fe x] profile of NGC 
4151 applies to NGC 5548 as well. The blue wing of the 
[Fe x] profile extends out to about 1500 km s - 1 relative to 
the systemic velocity (but about 1000 km s“ 1 relative to the 
peak). Even though v m . dX /v min may not be as large for this 


profile as other profiles presented thus far, a wide range of 
radial velocity is clearly present in this profile. Thus, the 
ionization parameter will be fairly constant over much of 
the NLR in NGC 5548. The blue wing of [O i] is about as 
extended as that of [Fe x] (i.e., significant flux out to a 
projected velocity of about — 1000 km s 1 ), from which we 
infer that the clouds are also optically thick throughout the 
NLR. 

Note that the [O i] profile is broad, but the [S ii] profile 
is narrow, in contrast to the general trend discussed here 
and in MC. As discussed further in § 5, for NGC 5548, the 
widths of [O i] 26300 and [O in] 24363, as well as the 
median velocity of both these lines, seem to correlate better 
with critical density than with ionization potential. Because 
NGC 5548 is unusual in other respects as well, we discuss 
this object in particular detail in § 6. 

5. PROFILE PARAMETER CORRELATIONS 

So far, the features of the profiles we have emphasized are 
usually obvious, and we have had little need for numerical 
values of profile parameters. For example, when we have 
noted that [O i] and [S ii] profiles are narrow compared 
with [Fe vn], the difference in width is large enough to 
discern immediately by visual inspection of the profiles. 
Indeed, much of the discussion of line width correlations 
that follows is driven by the similarity of [O i] and [S n] 
line widths as discussed in MC. However, profile param- 
eters are useful to confirm continuity in the sense that lines 
of intermediate critical density and/or ionization potential 
have intermediate profile parameters. 

We have adopted the area-defined parameters of Whittle 
(1982, 1985a) with a slight modification in the wings to 
allow for pixels that have negative flux after continuum 
subtraction. In Table 2 and Figures 8-13, we present corre- 
lations involving the Whittle parameters, critical density, 
and ionization potential. We have chosen the median veloc- 
ity, F(50%), and the interpercentile velocity width at 20%. 
The later is defined in Whittle (1982, 1985a) as IPV 
(20%) = a + 6, where a is the velocity (relative to the 
median) such that 10% of the total profile flux is blueward 
of a (see Fig. 7). Analogously, 10% of the total profile flux is 
redward of b . To be consistent with other studies, we have 
included analysis and discussion of emission-line blueshift. 
However, the sensitivity of the blueshift to details of the 
obscuration mechanism, which is poorly constrained, com- 
plicates the interpretation (see also Moore 1994, 1995). 
Therefore, we emphasize the line widths, because their cor- 
respondence to NLR kinematics is more direct than the 
blueshifts. 

For reference, we have also included in Table 2 the value 
of the IPV and asymmetry functions averaged between 20% 
and 80% (IPV avg and 4 avg ), the asymmetry function at 20% 
(A [20%]), and the kurtosis of the profiles. The asymmetry 
function is defined in Whittle ( 1 982, 1 985a) as 

«(%/ 2)-fo(%/2) 

' ° a(%/2) + b(%/2 ) ' 

The asymmetry function is positive for blueward extension 
and negative for redward extension. The kurtosis is defined 
in Whittle (1982, 1985a) as 2.5[IPV(50%)/IPV(10%)], where 
the normalization is chosen so that a Gaussian has kurtosis 
of 1. Kurtosis less than 1 corresponds to a profile more 
cusplike than a Gaussian, and the converse is true for 



TABLE 2 

Whittle Parameters 


Line 

K(50%) 
(km s “ l ) 

IPV(20%) IPV 

(kms M ) (km s _ ^) 

Kurtosis 

A(20%) 


MCG 8-11-11 

HP A4861 

124 

828 

352 

0.87 

0.09 

0.11 

[O hi] 25007 

92 

775 

340 

0.90 

0.07 

0.06 

[Fe vn] 25721 

16 

646 

301 

1.10 

0.15 

0.07 

He i 25876 

43 

665 

308 

1.05 

0.20 

0.11 

[Fe vn] 26087 

-29 

867 

371 

0.90 

0.18 

0.09 

[O i] 26300 

1 

871 

356 

0.80 

0.12 

0.05 

Ha 26563 

105 

454 

205 

1.03 

0.09 

0.05 

[S n] 26716 

78 

695 

299 

0.91 

0.08 

0.03 

Mrk 79 

[Fe Vn] 23760 

-141 

906 

427 

1.12 

0.05 

0.13 

[Ne hi] 23869 

53 

440 

190 

0.87 

0.14 

0.05 

HP 24861 

44 

400 

183 

0.98 

0.08 

0.06 

[O hi] 25007 

61 

493 

205 

0.69 

-0.05 

-0.01 

[Fe vn] 26087 

-30 

708 

295 

0.90 

0.30 

0.15 

[O i] 26300 

58 

403 

175 

0.84 

-0.04 

0.00 

[Fe x] 26374 

-113 

838 

360 

1.01 

0.26 

0.08 

Ha 26563 

35 

379 

170 

1.00 

0.00 

-0.01 

[N ii] 26584 

21 

424 

184 

0.93 

0.12 

0.03 

[S ii] 26731 

19 

313 

142 

1.06 

0.16 

0.08 

Mrk 704 

[Fe vn] 23760 

-122 

572 

231 

0.777 

0.232 

0.143 

[Ne in] 23869 

31 

449 

192 

0.866 

0.002 

-0.028 

[O hi] 25007 

62 

418 

183 

0.911 

0.061 

0.015 

[Fe vn] 25721 

-96 

412 

168 

0.834 

0.236 

0.114 

[Fe vn] 26087 

-94 

649 

254 

0.769 

0.291 

0.180 

[O i] 26300 

69 

294 

136 

1.13 

0.037 

0.018 

[Fe x] 26374 

-102 

479 

215 

0.96 

0.093 

0.040 

Ha 26563 

-1 

338 

154 

1.05 

-0.075 

-0.053 

[N ii] 26584 

80 

305 

138 

1.04 

0.030 

0.024 

[S ii] 26731 

52 

340 

154 

1.08 

0.146 

0.041 

Mrk 841 

[O hi] 24363 

78 

306 

128 

0.89 

-0.04 

-0.03 

He ii 24686 

70 

268 

114 

0.92 

-0.02 

-0.05 

HP 24861 

51 

421 

145 

0.69 

0.41 

0.15 

[O m] 25007 

66 

306 

119 

0.69 

0.06 

0.01 

He i 25876 

53 

165 

73 

1.02 

-0.02 

0.00 

[Fe vu] 26087 

-138 

1435 

543 

0.76 

0.44 

0.22 

[O I] 26300 

39 

193 

84 

0.94 

-0.05 

-0.02 

[Fe x] 26374 

-47 

394 

181 

1.08 

0.20 

0.17 

[N ii] 26548 

48 

139 

62 

1.02 

0.13 

0.06 

Ha 26563 

39 

306 

117 

0.72 

0.12 

0.04 

[N ii] 26584 

32 

280 

110 

0.78 

0.25 

0.10 

[S ii] 26716 

36 

186 

80 

0.92 

0.07 

0.02 

[S ii] 26731 

21 

201 

87 

0.89 

0.03 

-0.02 

NGC 4151 

[Fe vn] 23760 

-114 

721 

299 

0.86 

0.12 

0.09 

[Ne in] 23869 

-55 

111 

292 

0.71 

0.19 

0.11 

[Fe vn] 26087 

-49 

791 

285 

0.68 

0.26 

0.17 

[O i] 26300 

36 

541 

208 

0.74 

0.05 

0.04 

[Fe x] 26374 

-156 

1054 

500 

1.08 

0.29 

0.21 

Ha 26563 

32 

381 

155 

0.86 

0.03 

0.02 

NGC 5548 

[Fe vn] 23760 

-376 

564 

252 

1.03 

0.05 

0.04 

[Ne in] 23869 

-151 

574 

259 

0.97 

0.04 

0.03 

[O hi] 24363 

-303 

521 

225 

0.91 

0.24 

0.08 

HP 24861 

-98 

491 

224 

1.05 

0.08 

0.05 

[O m] 25007 

-132 

547 

234 

0.85 

0.14 

0.07 

[Fe vn] 25721 

-287 

572 

262 

0.99 

0.00 

0.02 

[Fe vn] 26087a 

-289 

717 

303 

0.83 

0.05 

0.02 

[Fe vn] 26087b 

-334 

726 

302 

0.89 

0.10 

0.02 

[O i] 26300 

-127 

627 

249 

0.76 

0.31 

0.16 

[Fe x] 26374 

-436 

762 

304 

0.82 

0.17 

0.00 

Ha 26563 

-57 

417 

195 

1.12 

0.12 

0.06 

[N ii] 26584 

31 

283 

126 

1.01 

0.13 

0.05 

[S ii] 26716 

26 

306 

144 

1.09 

0.07 

0.05 

[S ii] 26731 

19 

327 

149 

1.04 

0.07 

0.02 
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lb) 

Asymmetry function 


(c) 

IPV function 



Fig. 7. — Definitions for an area parameter scheme to characterize profile asymmetry and width, (a) 10, 50, and 90 interpercentile markers in the [O in] 
profile of NGC 5643 (wavelength increasing to the right). Lengths a and b are used to evaluate asymmetry and width parameters .4(20%) and IPV(20%), 
respectively. ( b , c) Asymmetry and IPV functions for this profile. This figure and caption are both reproduced from Whittle (1985a). We thank Mark Whittle, 
as well as the MNRAS, for permission to include them here. 






Fig. 8. — The median velocity, F(50%), and the interpercentile velocity width at 20% (see text and Fig. 7) are plotted vs. n cr and IP for MCG 8-11-11. 
Isolated open circles represent the recombination lines, and the open circle enclosing the filled circle represents [O i] A6300. There is no correlation of 
IPV(20%) with either n cr or IP. 
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Fig. 9. — The median velocity, F(50%), and the interpercentile velocity width at 20% (see text) are plotted vs. n CI and IP for Mrk 79. Isolated open circles 
represent the recombination lines. Although IPV(20%) correlates with both n cr and IP, [O i] 26300 (enclosed by an open circle) lies off the correlation on the 
n cr plot, but it follows the correlation on the IP plot. 


kurtosis greater than 1. The values of K(50%), IPV(20%), 
and IPV avg are in units of km s -1 , while the kurtosis, 
T(20%), and A avg are dimensionless. 

For possible future reference, we present six additional 
profile parameters in Table 3. The velocity of the peak and 
centroid are F peak and F ccn troid‘ The remaining parameters 
are widths as follows. The rms velocity is K rms = (<T 2 > 
— <u> 2 ) 1/2 . The full widths at and j of the peak inten- 
sity are FW^M, FW^M, and FWHM, respectively. All the 
parameters of Table 3 are in units of km s -1 . 

As an additional reference, we have also tabulated those 
emission lines we have analyzed together with their critical 
density and ionization potential in Table 4. As is conven- 
tional, we define the ionization potential as the arithmetic 
mean of the two successive stages of ionization. 

In Figures 8-13, we plot F(50%) and IPV(20%) in units 
of km s" 1 versus the log of n cr in units cm -3 and versus IP 
in eV for all six of our objects. Because much of the dis- 
cussion that follows focuses on [O i] 26300 and [O in] 
24363 as the lines that can potentially discriminate between 
correlations with n CT and correlations with IP, we highlight 


the data points corresponding to these lines by enclosing 
them with open circles. The isolated open circles represent 
Ha and Hj8. Although the critical density and ionization 
potential are not physically meaningful for the recombi- 
nation lines, the profiles contain information. In stratified 
clouds, the recombination lines are emitted most efficiently 
in a zone of intermediate ionization, in which [O m] and 
[Ne m] are also emitted. Due to the relatively large uncer- 
tainties of subtracting the underlying broad component 
from the narrow component, it is not appropriate here to 
attach much significance to the similar velocity shifts (and 
line widths) of the Balmer lines and the lines of [O hi] 25007 
and [Ne in] 23869. However, we suggest that the recombi- 
nation lines should be included in similar analyses of Seyfert 
2 galaxies. 

For NGC 5548, correlations of both F(50%) and 
IPV(20%) are clearly with n CI rather than with IP, making 
this object a counterexample to the general trend discussed 
here and in MC. Therefore, we discuss the correlations for 
NGC 5548 separately. For the other objects, correlations 
are less clear, particularly the correlation with K(50%). 









3 4 5 6 7 8 9 10 0 50 100 150 200 250 

Log of Critical Density (cm' 3 ) Ionization Potential (ev) 


Fig. 10. — The median velocity, F(50%), and the interpercentile velocity width at 20% (see text) are plotted vs. n CT and IP for Mrk 704. Isolated open circles 
represent the recombination lines. Although IPV(20%) correlates with both n CT and IP, [O i] 26300 (enclosed by an open circle) lies off the correlation on the 
n CT plot but follows the correlation on the IP plot. 


Ultimately, possible correlations between 7(50%) and/or not very sensitive to either the strength or the slope of the 

IPV(20%) with n cr and/or IP must be accepted or rejected synthetic continuum. The reason is that we know a priori 

by comparing visually the scatter in the result with some the functional form of the adjacent continuum, though this, 

systematic trend or lack thereof. Error bars that truly reflect of course, is not the case for the real data. The standard 

a quantitative uncertainty in the parameters in cases in deviations are also not very sensitive to the precise choice of 

which nonformal errors dominate, such as that of [Fe x] profile limits, though this again is at least partly an artifact 

26374, are virtually impossible to estimate reliably. Only in of the continuum subtraction being more reliable than it is 

the case in which the adjacent continuum is nearly flat or for the real data. 

linear, and featureless, and in which there are no blends Letting /be the noise added to the profile relative to the 

with other lines, are the true uncertainties approximated peak of the profile, the standard deviation in 7(50%) 

reasonably well by formal error estimates. increases linearly from 0 to 70 km s" 1 as/ is increased from 

We estimated formal errors by choosing the [O in] 25007 0 to 0.3. For IPV(20%), the standard deviation increases 

profile in NGC 5548 as a fiducial model profile. We added a linearly from 0 to 230 km s _1 . (The standard deviation 

linear continuum and noise to the profile, subtracted a curves are more complicated for/ > 0.3.) The noisiest pro- 

linear fit to the continuum, and calculated 7(50%) and files that we analyzed have/ < 0.1. It follows that the scatter 

IPV(20%) for the resulting profile. We repeated the pro- in the following plots that do not show correlations (as well 

cedure 10 times and calculated the standard deviation in the as those that do) is due mostly to either nonformal errors or 

10 calculated values of 7(50%) and IPV(20%). We repeated to the likelihood that the objects are more complicated than 

the experiment with varying amounts of noise and contin- justifies assuming an intrinsic one-to-one correspondence 

uum strength and slope. For fixed S/N in the final between line width and n CT and/or IP. In Paper II, for 

continuum-subtracted profile, the standard deviations are example, we present synthetic correlations resulting from 
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Fig. 1 1. — The median velocity, F(50%), and the interpercentile velocity width at 20% (see text) are plotted vs. n cr and IP for Mrk 841. Isolated open circles 
represent the recombination lines, and open circles enclosing the filled circles represent [O i] 26300 and [O hi] 24363. It is inconclusive whether correlations 
are with IP or n cr because only the [Fe] lines, which have both high n CI and high IP, stand out. 


models, in which scatter comparable to that observed is 
intrinsic to the model. 

Now we present and discuss the correlations for individ- 
ual objects. For each correlation, we also calculate a linear 
correlation coefficient, R , together with the corresponding 
probability of chance occurance, P null . 

MCG 8-11-11 (Fig. 8). — V(50%) appears to correlate 
better with n cr than with IP. However, because [O hi] >15007 
is not blueshifted relative to [S n] A6716, any correlation 
between blueshift and critical density that might be inferred 
from Figure 8a is at best tenuous. There is no correlation of 
IPV(20%) with either n CT or IP (we have already noted in 
§5.1 that the widths, as well as the profiles, of lines of both 
high and low IP are similar in this object). For MCG 8-11- 
11, R/P null are -0.73/0.16 for F(50%) versus n cr , -0.50/ 
0.39 for F(50%) versus IP, 0.27/0.66 for IPV(20%) versus 
u cr , and —0.18/0.78 for IPV(20%) versus IP, respectively. 

Mrk 79 (Fig. 9). — It appears that the highest critical 
density lines are blueshifted relative to the low critical 
density lines (Fig. 9a). Note, however, that these lines have 
high IP and that [O i] A6300 is associated with the lines of 


low-ionization potential rather than with those of high criti- 
cal density. Although a correlation with n cr is apparent 
visually, to interpret Figure 9a as a correlation of F(50%) 
with n cr , we must suppose that F(50%) is roughly constant 
over four decades in critical density and then changes 
rapidly over an additional factor of order unity in w cr . For 
IPV(20%), ^6300 lies off the correlation on the critical 
density plot (Fig. 9c) but follows the correlation on the 
ionization potential plot (Fig. 9 d). The values of R/P nuU are 
-0.64/0.08 for F(50%) versus n CTf -0.79/0.02 for F(50%) 
versus IP, 0.83/0.02 for IPV(20%) versus n CT , and 0.84/0.01 
for IPV(20%) versus IP, respectively. 

Mrk 704 (Fig. 10). — On the F(50%) correlation, the [O i] 
line fits on both the n cr and the IP relations. On the IPV 
plot, the [O i] line fits well on the IP correlation and falls 
out of place on the n cr correlation. The values of R/P nuU are 
-0.79/0.01, 0.61/0.08, -0.82/0.006, and 0.58/0.10. 

Mrk 841 (Fig. 11). — For all the correlations for Mrk 841, 
only the [Fe vii] and [Fe x] lines stand out. It is inconclu- 
sive whether correlations are with IP or with w cr because 
these lines have both high IP and high n cr . The values of 
R/P null are -0.46/0.22, -0.62/0.07, 0.44/0.24, and 0.40/0.28. 







TABLE 3 

Additional Profile Parameters 


Line 


y 

r centroid 

v 

r rtm 

FWjM 

FW^M 

FWHM 

MCG 8-11-11 

[O ii] 23727 

71 

-46 

395 

1134 

1011 

823 

H/9 >1486 1 

115 

48 

331 

888 

764 

544 

[O m] 25007 

119 

19 

323 

899 

757 

624 

[Fe vn] 25721 

126 

-50 

279 

895 

811 

607 

He i 25876 

83 

-26 

288 

810 

748 

650 

[Fe vii] 26087 

136 

-106 

351 

955 

846 

606 

[O i] 26300 

117 

-74 

390 

912 

755 

577 

Ha 26563 

121 

56 

197 

611 

544 

413 

[S ii] 26716 

64 

30 

287 

1185 

690 

555 

Mrk 79 

[O ii] 23727 

23 

-46 

201 

725 

658 

550 

[Fe vii] 23760 

-160 

-249 

318 

1009 

1003 

854 

[Ne hi] 23869 

-12 

-49 

179 

528 

464 

433 

H/J 24861 

18 

-30 

159 

509 

451 

367 

[O hi] 25007 

-7 

25 

260 

506 

452 

425 

[Fe vii] 26087 

-4 

-119 

316 

733 

593 

415 

[O i] 26300 

25 

53 

241 

456 

400 

353 

[Fe x] 26374 

-31 

-201 

306 

906 

772 

570 

Ha 26563 

-8 

3 

171 

458 

422 

342 

[N ii] 26584 

-16 

-27 

195 

471 

425 

344 

[S ii] 26731 

13 

-32 

132 

409 

383 

313 

Mrk 704 

[O ii] 23727 

92 

43 

179 

737 

597 

448 

[Fe vii] 23760 

-175 

-257 

190 

627 

402 

298 

[Ne m] 23869 

-70 

-50 

183 

512 

446 

362 

[O iii] 25007 

0 

-12 

161 

476 

429 

393 

[Fe vii] 25721 

-103 

-158 

185 

483 

315 

233 

[Fe vii] 26087 

-94 

-181 

303 

511 

410 

262 

[O i] 26300 

42 

30 

256 

416 

379 

315 

[Fe x] 26374 

-116 

-155 

217 

544 

576 

460 

Ha 26563 

-67 

-29 

148 

457 

401 

318 

[N ii] 26584 

60 

42 

135 

405 

365 

292 

[S ii] 26731 

104 

26 

137 

473 

439 

305 

Mrk 841 

[O ii] 23727 

-55 

31 

203 

552 

515 

445 

[O m] 24363 

-3 

7 

112 

394 

347 

218 

He ii 24686 

-12 

13 

102 

351 

280 

206 

HP 24861 

4 

-53 

155 

300 

283 

212 

[O hi] 25007 

0 

-2 

153 

296 

258 

189 

He i 25876 

14 

20 

71 

246 

197 

148 

[Fe vii] 26087 

3 

-351 

588 

1337 

718 

639 

[O i] 26300 

0 

6 

90 

233 

204 

166 

[Fe x] 26374 

-58 

-112 

160 

491 

438 

313 

[N ii] 26548 

21 

12 

61 

197 

169 

128 

Ha 26563 

9 

-5 

148 

267 

231 

171 

[N ii] 26584 

10 

-20 

127 

258 

216 

164 

[S ii] 26716 

4 

-2 

87 

224 

192 

150 

[S ii] 26731 

1 

4 

89 

224 

203 

167 

NGC 4151 

[O ii] 23727 

-137 

-116 

301 

795 

723 

543 

[Fe vii] 23760 

-105 

-217 

280 

988 

907 

678 

[Ne iii] 23869 

-107 

-107 

320 

644 

511 

370 

[Fe vii] 26087 

-2 

-89 

372 

520 

390 

283 

[O i] 26300 

4 

-33 

231 

500 

391 

245 

Fe x 26374 

-2 

-233 

437 

1092 

995 

543 

Ha 26563 

-1 

-7 

166 

443 

353 

240 

NGC 548 

[Fe vii] 23760 

-430 

-430 

224 

1478 

1141 

490 

[Ne iii] 23869 

-149 

-243 

236 

682 

610 

507 

[O in] 24363 

-325 

-387 

224 

677 

573 

470 

Hp 24861 

-129 

-174 

182 

621 

588 

484 

[O hi] 25007 

-164 

-227 

248 

605 

538 

411 

[Fe vii] 26087 

-325 

-323 

324 

790 

675 

545 

[O i] 26300 

-120 

-219 

275 

542 

464 

329 

[Fe x] 26374 

-556 

-494 

301 

777 

633 

500 

Ha 26563 

-65 

-105 

179 

560 

517 

424 

[N ii] 26584 

33 

-12 

121 

378 

323 

243 

[S ii] 26716 

17 

-11 

127 

407 

366 

300 

[S n] 26731 

-19 

-21 

131 

422 

393 

307 
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TABLE 4 


Emission Lines, 2 , n cr , and IP 


Wavelength 

(A) 

Ion 

n cr 

(cm -3 ) 

Mean Ionization 
Potential 
(eV) 

3726.1 

[On] 

3.2 x 10 3 

24 

3728.8 

[On] 

6.3 x 10 2 

24 

3760.3 

TFe vii] 

4.0 x 10 7 

114 

3868.8 

fNe ml 

1.0 x 10 7 

52 

4363.2 

[Om] 

3.2 x 10 7 

45 

4685.7 

He ii 


40 

4861.3 

H i 



5006.8 

[O in] 

6.3 x 10 5 

45 

5721.1 

[Fe vii] 

4.0 x 10 7 

114 

5875.7 

He i 


12 

6086.9 

[Fe vii] 

4.0 x 10 7 

114 

6300.3 

[O i] 

2.0 x 10 6 

7 

6374.6 

[Fe x] 

5.0 x 10 9 

248 

6548.1 

[Nil] 

7.9 x 10 4 

22 

6562.8 

H i 



6583.4 

[Nil] 

7.9 x 10 4 

22 

6716.4 

[S ii] 

1.6 x 10 3 

17 

6730.8 

[S II] 

4.0 x 10 3 

17 



NGC 4151 (Fig. 12). — Because we do not have [S n] 
26716, 6731 profiles from NGC 4151 (and [O ii] 23727 is an 
unresolved blend), all we can say about this object, based on 
our data alone, is that the highest IP/n cr lines are blue- 
shifted and broad relative to [O i] 26300. However, if we 
were to include the [S ii] profiles of Veilleux (1991a), we 
would put NGC 4151 in the same category as MCG 8-11- 
1 1 (no correlation of IPV with either ionization potential or 
critical density). 

NGC 5548 (Fig. 13). — Because the [O i] profile in NGC 
5548 is broad and that of [S ii] is narrow, NGC 5548 stands 
in striking contrast to the general trend discussed here and 
in MC. In Figures 13a and 136, we plot the median velocity 
(F[50%]) of the emission lines from NGC 5548 versus both 
the critical density and the ionization potential, respec- 
tively. The velocities are relative to a 21 cm measurement of 
the systemic velocity. 

As stated already, for NGC 5548, correlations of both 
K(50%) and IPV(20%) are clearly with n CT rather than with 
IP. There is a good correlation between blueshift and criti- 





Fig. 12. — The median velocity, K(50%), and the interpercentile velocity width at 20% (see text) are plotted vs. n cr and IP for NGC 4151. Isolated open 
circle represents Ha, and open circle enclosing the filled circle represents [O i] 26300. If we were to include the [S n] 226716, 6731 lines of Veilleux (1991a), we 
would infer no correlation of IPV(20%) with either n CT or IP. 
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Fig. 13. — The median velocity, K(50%), and the interpercentile velocity width at 20% (see text) are plotted vs. n cr and IP for NGC 5548. Isolated open 
circles represent the recombination lines. The important [O i] A6300 and [O hi] A4363 data points are enclosed by open circles. Correlations of both F(50%) 
and IPV(20%) are clearly with n CT rather than with IP. This is in contrast to the general trend discussed here and in MC. 


cal density in NGC 5548. Although there is apparently a 
good correlation between blueshift and ionization potential 
as well, there are two points on the K(50%) versus ioniza- 
tion potential plot that deviate. That these deviant points 
represent [O i] 26300 and [O hi] 24363 is not surprising 
because, of the lines commonly analyzed in the literature, it 
is these two that deviate from a correlation of critical 
density with ionization potential. 

Note that the velocity shift from systemic is very large 
(^400 kms -1 ) for the forbidden Fe lines in NGC 5548. The 
shift is large enough that one may suspect an error in the 
wavelength calibration. This is not the case, however, 
because the wavelength calibration matches the known 
wavelengths of the night-sky lines to a precision of 25 km 
s -1 . One might also speculate that the 21 cm systemic 
velocity measurement is not representative of the NLR 
center-of mass velocity. However, the 21 cm systemic veloc- 
ity agrees (within 40 km s 1 ) with F(50%) of the lowest 
critical density lines. Regardless of the systemic velocity, the 
Fe lines are blueshifted by about 400 km s -1 relative to the 


lowest critical density lines. Thus, NGC 5548 is also 
unusual in its large blueshift in F(50%) of the Fe lines. 

In Figures 13c and 13 d, we plot IPV (20%) versus, n cr and 
IP. Again, [O i] 26300 deviates from the otherwise well- 
defined correlation on the IP plot but does not deviate on 
the n cr plot. The above correlations for NGC 5548 are con- 
sistent with velocity decreasing outward because lines of 
high critical density are systematically broader than lines of 
low critical density. In our previous discussion, we have 
assumed that all Seyfert NLRs are fundamentally similar 
and can be described by models differing only in, e.g., 
N c JN min at v min . Although the line profiles (especially that 
of [Fe x] 26374) are similar to the profiles of the other 
objects in our sample, the correlations suggest an NLR with 
fundamentally different kinematics. In MC, we suggested 
that at 5 A and 10 A resolution, 24363 line widths may have 
been systematically overestimated in Seyfert 1 galaxies 
because of difficulties in deblending from the broad and 
narrow components of Hy. Even with our higher resolution, 
we may be subject to the same effect as well. However, the 
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dilemma in which line widths of NGC 5548 seem to correl- 
ate fundamentally with critical density rather than with ion- 
ization parameter is posed by both A4363 and 26300. The 
(highly significant) values of R/P null are —0.94/5.0 x 10” 6 , 
-0.83/9.2 x 10” 4 , 0.88/1.8 x 10“ 4 , and 0.71/9.7 x 10" 3 . 

6. SUMMARY 

The most important observational feature of the profiles 
presented here is that for lines encompassing a wide range 
of critical density and ionization potential, the variation of 
flux with projected velocity is very similar in some objects 
(e.g., MCG 8-11-11 and NGC 4151). This is inconsistent 
with a scenario in which lines of low and high critical 
density and/or ionization potential are formed in physically 
distinct regions. Thus, we expect comparable spatial extents 
for emission into each line. We have proposed testing this 
hypothesis with coordinated Hubble Space Telescope and 
ground-based spectrophotometry of the applicable NLRs 
that can be spatially resolved on an angular scale ( < 2") that 
properly complements analyses of existing spatially inte- 
grated emission-line spectra. In some other objects, lines of 
high critical density and/or ionization potential are system- 
atically broader than those of low critical density and/or 
ionization potential (e.g., Mrk 79, Mrk 704, and Mrk 841). 
It is possible to explain the difference between the two types 
of objects in terms of a relatively small difference in the 
cloud column density throughout the NLR (see Paper II, 
MC, and Moore [1994]). If clouds at small r have N col % 
10 24 cm -2 , then even with N C0l cc r~ 413 (constant mass 
spherical clouds with n e oc r~ 2 ), it is possible to have N co , > 
N mm ( ~ 10 22 cm -2 for T = 0.3) throughout the NLR so 
that each cloud emits both [O i] and [Fe vii]. If clouds start 
out with Af col less by a modest factor, then iV co , falls to less 
than N m in before clouds reach their maximum velocity at 
large r, and an [O i] profile narrower than that of [Fe vii] 
results. 


Overall, our sample of Seyfert 1 galaxies confirms the 
well-known correlations of line width and blueshift with 
ionization potential and/or critical density. A sample the 
size of ours cannot alone confirm or refute properly the 
conclusion of MC that the correlations are fundamentally 
with ionization potential rather than with critical density, 
especially given that our sample includes one unusual 
counter example (NGC 5548). The correlations for Mrk 841 
are not conclusive because the only lines that stand out are 
[Fe vii] 26087 and [Fe x] >16374, which have both high 
critical density and high ionization potential. This leaves 
just four objects, MCG 8-11-11, Mrk 79, Mrk 704, and 
NGC 4151, to support our earlier conclusion (correlations 
with ionization potential). Nevertheless, these four objects 
do follow the trend discussed in MC; the objects with broad 
[O i] profiles (MCG 8-11-11 and NGC 4151) also have 
broad [S ii] profiles, and the objects with narrow [O i] 
profiles (Mrk 79 and Mrk 704) have narrow [S n] profiles. 
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